Internal stucture of clusters in 112 122 Ba nuclei within relativistic mean field theory 
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We study the clustering structure and the internal or sub- structure of clusters in 112-1 22 Ba nuclei within the 
framework of relativistic mean field theory in an axially deformed cylindrical co-ordinate. We calculate the 
total density distribution, and the individual neutrons and protons density distributions. From the analysis of the 
clustering confugurations of the density distributions of various shapes, we find different sub- structures inside 
the Ba nuclei considered here. The important step, carried out for the first time, is the counting of number of 
protons and neutrons present in the clustering region(s). 12 C is shown to consitute the cluster configuration of 
Ba nuclei in most cases, with 2 ' 3 H and 4 He constituting the neck between two fissioning symmetrical fragments. 

PACS numbers: 21.10.Dr, 21.60.-n, 21.60.Jz, 21.60.Gx, 24.10.Jv, 23.70.+j 



I. INTRODUCTION 

Instabilities of "stable" nuclei against exotic cluster de- 
cays were first studied by one of us (RKG) and collaborators 
Q S H for the neutron-deficient 108 - 116 Xe, 112 - 120 Ba, 

116-124 Ce9 120-124 Nd? 124-128 Sm? ^~^Qd and 

neutron-rich 146 Ba, 152 Ce, 156 Nd, 160 Sm, and 164 Gd nu- 
clei, on the basis of the preformed cluster model (PCM) of 
Malik and Gupta [H, 0] (see, also the works of other au- 
thors 0, Si). The N = Z, a-nuclei 8 Be, 12 C, 16 0, 20 Ne, 
24 Mg, and 28 Si are predicted to be the most probable emit- 
ters from neutron-deficient parents, with a much smaller de- 
cay half-life compared to N > Z clusters predicted prefer- 
ably or observed from neutron-rich parents (e.g., 14 C from 
146 Ba or 222 Ra with 132 Sn or 208 Pb as daughter, respectively). 
Thus, other than a particle, 12 C decay of 112_120 Ba nuclei is 
shown to be the most probable one with 100 Sn and its heav- 
ier isotopes as the daughter nuclei. The ground- state decay 
of Ba, however, could not be established as yet |H, 0U, and 
a new phenomenon of intermediate mass fragments (IMFs) 
with 3<Z<9, also referred to as 'clusters' or 'complex frag- 
ments', is observed to be emitted from the compound nu- 
cleus 116 Ba* formed in 58 Ni+ 58 Ni^ 116 Ba* reactions at both 
the high [Qjl [13] and medium energies (TH [3]. The mea- 
sured IMF cross-section cfimf for the 116 Ba* decay at all the 
above mentioned medium and high energies are so far under- 
stood only on the preformed-clusters (PCM) based dynamical 
cluster-decay model (DCM) of one of us (RKG) and collab- 
orators lfl5ll . The DCM describes the 116 Ba* data on cf IM f 
reasonably well, and predicts an additional fusion-fission of 
116 Ba* which consists of fragments at the heavy end of sym- 
metric and near symmetric division (14<Z<28), very recently 
observed at GANIL [[H] for the decays of 118 > 122 Ba* nuclei 
produced in 78 ' 82 Kr+ 40 Ca reactions at a lower incident en- 
ergy, and also explained more recently B17I1 on DCM. It should 
thus be interesting to see if mean-field approaches, such as the 
relativistic mean field model, support such a clustering struc- 
ture for the ground and/ or excited states of Ba nuclei. 

Relativistic mean field (RMF) calculations for light nuclei 
frill have very successfully shown the a — a clustering in 8 Be, 



a benchmark nucleus, and the a-clustering and halo structures 
of (a+a)-core or (a+a+p)-core plus xn for the stable and ex- 
otic Be and B nuclei. For the N=Z, a-nuclei 12 C and 16 O, the 
RMF formalism gives the several known ground- state (g.s.) 
structures such as the 3a-equilateral triangle (co-existing with 
spherical shape) and 4a-tetrahedron or kite-like, plus the 3a- 
and 4a-linear chains for their excited states. Also, both a- and 
non-a-cluster structures, like 10 B+ 10 B, 12 C+ 12 C, 12 C+a+ 12 C 
and 16 0+ 16 0, and 5-a-trigonal bipyramid and pentagonal 
bipyramid (hollow sphere) are obtained for 20 Ne to 32 S nu- 
clei, but no 5-a, 6-a chains, etc., are seen. Furthermore, 56 Ni 
shows lll9ll a preferred N=Z, a-nucleus clustering for states 
with deformations up to hyper-deformation ((3 <2.45). Simi- 
larly, for heavy actinides ( 222 Ra, 232 U, 236 Pu and 242 Cm), the 
RMF gives [ 20] not only the N~ Z, a-like clustering in the 
g.s. but also the exotic N^Z clustering in excited states. Sig- 
natures of clustering structure in RMF calculations of super- 
heavy Z=114 and 120, N=172-184 nuclei (U [H are also 
obtained in terms of exotic N^Z clusters at the centre of their 
super-deformed g.s. or the clustering in to two large and some 
small pieces is universal for all super-deformed ground states 
in Z=120 nuclei. The super- superheavy 238 TJ+ 238 U^ 476 184* 
nucleus also supports the clustering phenomenon, with a kind 
of triple fission of an exotic cluster in the neck region of two 
equal fragments of N=Z matter II 2 211 . The actual internal or 
sub- structure structure of the clusters, however, was not deter- 
mined in these calculations, which is one of the aims of the 
present study. 



Taking RMF model as an established tool for the cluster 
structure in nuclei, in this paper we demonstrate its applica- 
tion to medium-heavy 112_122 Ba nuclei, with an additional 
attempt to determine for the first time the internal structure(s) 
of cluster(s), i.e., the number of protons and neutrons in a clus- 
ter. The paper is so designed that Section II describes the rel- 
ativistic mean field theory, and Section III gives the results 
obtained from our calculation. Finally a brief summary and 
concluding remarks are given in Section IV. 
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II. THE RELATIVISTIC MEAN FIELD (RMF) METHOD 

The relativistic Lagrangian density for a nucleon-meson 
many -body system l!23L l24ll 

C = Ji{i^3^ - Mjijji + ^ad„a - \mlcj 2 

+\m 2 w v*v li + \c 3 (v^vn 2 - gwQa^iV^ 

--B^.B^ + ]^m 2 p W.R^ - grfa^i-W* 
—F^F^ - e^ Q—p^ frAp. (1) 

From this Lagrangian we obtain the field equations for the 
nucleons and mesons. These equations are solved by ex- 
panding the upper and lower components of the Dirac spinors 
and the boson fields in an axially deformed harmonic oscil- 
lator basis with an initial deformation /3 . The set of cou- 
pled equations is solved numerically by a self-consistent it- 
eration method. The centre-of-mass motion (cm.) energy 
correction is estimated by the usual harmonic oscillator for- 
mula E c , mm = | (41 A -1 / 3 ). The quadrupole deformation pa- 
rameter /?2 is evaluated from the resulting proton and neutron 

quadrupole moments, as Q = Q n + Q p = yj AR 2 fa) . 

The root mean square (rms) matter radius is defined as (r^) = 
\ J p(r±, z)r 2 dr; here A is the mass number, and p(r±, z) 
is the deformed density. The total binding energy and other 
observables are also obtained by using the standard relations, 
given in I124I1 . We use the well known NL3 parameter set Il25n . 
This set reproduces the properties of stable nuclei. As out- 
puts, we obtain different potentials, densities, single-particle 
energy levels, radii, quadrupole deformations and the binding 
energies. For a given nucleus, the maximum binding energy 
corresponds to the ground state and other solutions are ob- 
tained as various excited intrinsic states. 

The density distribution of nucleons play the prominent role 
for studying the internal structure of a nucleus. For a differ- 
ent quadrupole deformation, the density distribution p(r±, z) 
inside the nucleus must vary. For example, the p(r±, z) for a 
spherical nucleus is symmetrical in (p, z)— plane. However, it 
is highly asymmetric for a largely deformed nucleus. Know- 
ing the density distribution of the spherical or (oblate/ prolate) 
deformed confuguration, we can calculate the number of nu- 
cleons for each confuguration, defined as 

n= / p(r±,z)dr, (2) 

J z\ J r\ 

with n as number of neutrons N, protons Z or mass A. Though 
a straight forward calculation, this is being carried out here for 
the first time. 



III. CALCULATIONS AND RESULTS 

A. Ground state properties of Ba nuclei using the RMF 
formalism 

1. Potential Energy Surface 

We first calculate the potential energy surface (PES) by us- 
ing the RMF formalism in a constrained calculation Jlfl H2I 
[2SL [30], i.e., instead of minimizing the Hq, we have min- 
imized H' = Ho — XQ2, with A as a Lagrange multiplier 
and Q2, the quadrupole moment. The term Ho is the Dirac 
mean field Hamiltonian for RMF model (the notations are 
standard and its form can be seen in Refs. JM M HI). 
In other words, we get the constrained binding energy from 

BE c = E« ^'S^ 1 ^ and the free ener gy from 
BE = J2ij < <^ J |^^> > • ^e converged free energy solu- 
tion does not depend on the initial guess value of the basis 
deformation (3o as long as it is nearer to the minimum in PES. 
However, it converges to some other local minimum when (3q 
is drastically different, and in this way we evaluate a different 
intrinsic isomeric state for a given nucleus. 
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FIG. 1: The potential energy surfaces for 112 Ba and 122 Ba nuclei, 
i.e, constrained binding energy (BE C ) as a function of quadrupole 
deformation parameter in RMF(NL3) calculations. 

The PES for the representative 112 Ba and 122 Ba nuclei are 
shown in Fig. 1 for a wide range of P2 starting from oblate 
to prolate deformation. The upper pannel is for 112 Ba and 
the lower one for 122 Ba. For 112 Ba, we notice the min- 
ima are around /? 2 =-0.39, -0.20, 0.23, 1.02 and 1.20, corre- 
sponding to binding energy BE C = 890.6, 894.2, 896.7, 884.2 
and 882.7 MeV, respectively, and their energy differences be- 
tween the nearest consecutive minima are 3.60, 2.50, 12.93 
and 1.52 MeV. However, in case of 122 Ba, only two minima 
exist around /3 2 =-0.23 and 0.26 at BE C = 1013.4 and 1012.4 
MeV. The intrinsic energy difference between these two con- 
figurations is 1.03 MeV. From this figure, it is clear that there 
exists two identical major minima at ~-0.23 and 0.23 in 
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TABLE I: The RMF(NL3) results of binding energy, the quadrupole deformation parameter 02, and the root mean scmare radii of charge r c 
and matter r m for the 112_122 Ba nuclei. The energy is in MeV, and radii in fm. The experimental data are from I3ll. 13211. and the binding 
energies marked with star (*) are the extrapolated values. The finite range droplet model (FRDM) binding energy and 02 are taken from 11331] 
and the theoretical charge radius r c is from Hartree-Fock-BCS (HFBCS) method [34]. No experimental data are available for r c . 



Nucleus 


BE 


Trms • 


RMF(NL3) 


P2 


BE 


r c 


02 




RMF(NL3) 


Expt. 




r c 


RMF(NL3) 


Expt. 


FRDM 


HFBCS 


FRDM 


A1 *Ba 


860.79 
893.50 




4.80 
4.65 


4.99 
4.79 


-0.39 
-0.17 












895.38 




4.62 


4.74 


0.24 




894.89 


4.72 


0.21 




882.36 




5.35 


5.48 


1.24 












897.33 




12.17 


12.21 


10.71 










114 Ba 


918.10 




4.78 


4.90 


-0.39 












920.13 


922.26 


4.65 


4.75 


0.24 




921.27 


4.74 


0.24 




909.35 




5.33 


5.45 


1.19 












920.61 




12.16 


12.20 


10.58 










116 Ba 


943.66 


947.024* 


4.80 


4.90 


-0.39 












945.70 




4.71 


4.79 


0.36 




946.85 


4.78 


0.28 




934.15 




5.36 


5.46 


1.20 












942.38 




12.13 


12.16 


10.40 










118 Ba 


969.16 




4.73 


4.81 


-0.24 












971.38 


971.022* 


4.75 


4.82 


0.33 




970.75 


4.80 


0.29 




962.45 




12.11 


12.13 


10.24 










120 Ba 


991.64 




4.75 


4.81 


-0.23 












993.94 


993.600 


4.77 


4.83 


0.32 




993.44 


4.82 


0.28 




975.04 




12.55 


12.58 


10.97 










122 Ba 


1013.53 




4.77 


4.82 


-0.22 












1015.52 


1015.499 


4.80 


4.84 


0.32 


0.354 


1015.21 


4.81 


0.27 




999.40 




12.09 


12.10 


9.95 











both the 112 Ba and 122 Ba nuclei. A further investigation of 
the diagram shows that actually the multi-minima structure of 
112 Ba disappears gradually with the increase of mass number 
in the isotopic chain of Ba, and reaches to only two minima 
configuration, one oblate and another prolate, at mass number 
A=122. 



2. Binding Energies and Quadrupole Deformation Parameter 

The binding energies and quadrupole deformation param- 
eters 02 for 112-1 22 Ba isotopes are evaluated for the ground 
as well as intrinsic excited states. The obtained results are 
tabulated in Table I, together with the experimental or extrap- 
olated values, wherever available. The experimental value of 
g.s. 02 is available for 122 Ba only. Since the neutron-deficient 
112-1 20 g a i so t p es lie ne ar the proton drip-line, their defor- 
mation parameters are not yet known. Table I shows that both 
the g.s. binding energies, charge radius (r c ) and 02 values 
agree well with the experimental data and with the theoretical 
calculations [133ll34j1. 

As discussed above, lighter Ba isotopes have several intrin- 
sic minima, where each minimum corresponds to a deforma- 
tion and a binding energy. The largest binding energy min- 
imum corresponds to the ground state and all other minima 
are the excited intrinsic states. In this way, the ground state 
02 is 0.24 for 112 Ba. Similarly, the g.s. deformations for 
114 Ba, 116 Ba, 118 Ba, 120 Ba and 122 Ba are 0.24, 0.36, 0.33, 
0.32 and 0.32, respectively. All other intrinsic state deforma- 



tions are also listed in Table I. For 112 ' 114 ' 116 Ba nuclei we 
get a solution at a highly deformed configuration of 02 ~ 1.2, 
whereas this hyper-deformed minimum is washed out with in- 
crease of mass number in the Ba isotopic chain. This means 
there is no hyper-deformed solutions for 118_122 Ba (see also, 
Fig. 4 where such highly deformed configurations are shown 
only for 112 ' 114 ' 116 Ba). If the nucleus is further deformed 
(02 ~10), it breaks into two fragments (see also Fig. 5) 
with a separation radius r m of about 12 fm. An analysis of 
this structure shows two distinct (fission-like) fragments con- 
nected with a thin density distribution of nucleons, which can 
be considered like a very thin neck. We find that the compo- 
sition of the (connecting) neck region is an isotopic chain of 
hydrogen or helium nuclei 2 ' 3 H and 4 He (discussed below). 



3. The root-mean-square radii r rms 

In this subsection, we discuss our calculation on root 
mean square (rms) radii of matter distribution, with various 
quadrupole deformations from oblate to prolate and then from 
super-deformation to hyper-deformation. The matter rms ra- 
dius r m and charge distribution rms radius r c are listed in Ta- 
ble I for different 02 values. From Table I, it is observed that 
the r m increases with increase of quadrupole deformation. Fi- 
nally, both r m and r c show a large extension of about 12 fm 
at 02 ~10.5. As already pointed out above, at this point of 
deformation, the two fragments get separated from each other 
leaving a thin density distribution of matter in between the two 
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fragments or clusters. 
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FIG. 2: (Color on line) The cluster confugurations of 112 122 B a for 
the oblate states in RMF(NL3). 
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FIG. 3: (Color on line) The cluster confugurations of 112 122 B a for 
the ground states in RMF(NL3). 



B. Cluster confugurations in 
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FIG. 4: (Color on line) The cluster confugurations of 112 122 Ba for 
the prolate states in RMF(NL3). 
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Figures 2 to 5 show the density plots for all the possible 
solutions, starting from oblate to prolate deformations of the 



FIG. 5: (Color on line) The cluster confugurations of 112 122 B a for 
the hyper-deformed states in RMF(NL3). 



Ba isotopic chain. Fig. 2 gives the density distributions for 
H2-i22g a f Qr ^ bl a te state solutions around the deforma- 
tion parameter ~-0.35, and Figs. 3, 4 and 5 show, re- 
spectively, the prolate state solutions of 112_122 Ba around 
/?2 ~0.25, 1.1 and 10.5. From these graphs the clustering 
structures in these nuclei are quite evident. Considering the 
colour code, starting from deep black with maximum to white 
(or light gray) bearing the minimum density, one can analyse 
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the distribution of nucleons inside the various isotopes at vari- 
ous shapes. In Fig. 2 the minimum density for the oblate- state 
starts from 0.001 fm -3 and goes up to maximum 0.16 fm -3 , 
but in case of the prolate- states in Figs. 3, 4 and 5, it starts 
from 0.001 fm -3 and goes up to 0.18 fm -3 , which shows that 
the size of cluster nucleus is smaller in the oblate- state than 
that of the prolate- state. A careful inspection of the density 
distributions in different regions of the nucleus clearly shows 
the formation of various cluster(s) inside the nuclei, which are 
listed in Table II as clusters I, II or III for a given deformation 

Further, if we increase the value of deformation parameter 
fa to ~ 10.4- 10.7 or more, the nucleus gets separated into two 
fragments with a neck formation as shown in Fig. 5 for the 
isotopic chain of 112_122 Ba nuclei. Another prominent obser- 
vation is that there is no configuration of deformation param- 
eter fa ~1.2 in case of 118_122 Ba, similar to what is shown in 
Fig. 4 for lighter mass Ba isotopes. If we compare the cluster 
confugurations of each nucleus, it is clear that there is a grad- 
ual change in the confuguration inside the nucleus, i.e., the 
clusters inside the nucleus for different fa values are distinct 
from each other, and for each isotope of 112_122 Ba chain. 



TABLE II: The number of nucleons A c i us ., the protons Z c i us ., and 
neutrons N c i us . in a cluster inside the 112-1 22 Ba nuclei for different 
solutions of deformations 02 obtained from the RMF(NL3) formal- 
ism. The range of the cluster (n, r<i \ zi, z?) are in fm. 



Nucleus 


02 


Cluster No. 
range (n, ri\z\, 22) 


A-clus. 


Zclus. 


Ndus. 


Cluster 


li2 Ba 


-0.39 


I(-0.7,0.7;-3.0,-1.7) 


11.6 


5.6 


6.0 


i2 C 




0.24 


I(1.9,4.5;-1.5,1.5 


36.0 


17.7 


18.3 


36 Ar 






II (-1.3,1.3;-2.6,-1.2) 


13.0 


6.3 


6.7 


13 C 




10.71 


I (-6.3,6.3;-3.5,3.5) 


1.6 


0.7 


0.8 


2 H 


114 Ba 


0.24 


I(2.1,4.6;-1.5,1.5) 


34.7 


16.9 


17.9 


35q 






II (-1.1,1.1;1.2,2.5) 


12.6 


6.1 


6.4 


13 C 






III (0.7,2.0;3.3,3.7) 


2.3 


1.1 


1.2 


2 H 




1.19 


I (3.9,5.4;-0.7,0.7) 


8.4 


4.1 


4.3 


8 Be 




10.58 


I (-6.6,6.3;-3.8,3.8) 


2.2 


1.1 


1.1 


2 H 


116 Ba 


0.36 


I(4.4,6.2;-1.6,1.6) 


25.3 


12.3 


13.1 


25 Mg 






II (2.6,3.4;-0.9,0.9) 


6.1 


2.8 


3.4 


6 Li 




10.40 


I(11.0,12.2;-0.8,0.8) 


7.5 


3.7 


3.8 


8 Be 






II (9.9,12.8;2.3,2.9) 


6.6 


3.2 


3.4 


6 Li 






III (-6.1,6.1;-3.4,3.4) 


1.9 


0.9 


1.0 


2 H 


118 Ba 


-0.24 


I(-0.6,0.6;-1.2,1.2) 


12.85 


6.9 


5.9 


13 C 




0.33 


I(2.5,5.5;-1.6,1.6) 


42.0 


19.9 


22.2 


42 Ca 






II (0.7,2.0;3.3,3.6) 


1.7 


0.8 


0.9 


2 H 






III (-0.4,0.4;-2. 1,-1.6) 


1.6 


0.7 


0.8 


2 H 




10.24 


I (-6.2,6.2;-3.7,3.7) 


2.8 


1.2 


1.6 


3 H 


120 Ba 


-0.23 


I(-0.6,0.6;-1.2,1.2) 


12.9 


5.8 


7.1 


13 C 




0.32 


I(-5.5,-5.0;-0.7,-1.2) 


0.9 


0.3 


0.6 


X H 




10.97 


I (-6.3,6.3;-4.3,4.3) 


3.8 


1.6 


2.1 


4 He 


122 Ba 


-0.22 


I(-1.0,1.0;-1.5,1.5) 


24.4 


11.2 


13.2 


23 Na 




0.32 


I(2.5,5.3;-1.7,1.7) 


42.6 


19.7 


22.8 


43 Ca 






II (-2.0,-0.8;3.4,3.7) 


1.6 


0.7 


0.8 


2 H 




9.96 


I-6.1,6.1;-3.7,3.7) 


2.4 


1.0 


1.4 


2 H 



C. Counting of nucleons in clusters formed inside the 
112122 Ba nuclei 



In this subsection, we count the number of nucleons in dif- 
ferent clusters formed inside the 112-122 Ba isotopic chain, 
listed in the Table II. The density distributions of these clus- 
ters, obtained from the RMF(NL3) formalism for different so- 
lutions of deformation parameters fa from oblate to prolate 
confugurations, are already shown above in Figs. 2 to 5. From 
these distribution plots of each oblate or prolate configuration, 
we find the number of nucleons by using the general formula 
in Eq. (2), for not only the total nucleons (using total density 
distribution of the nucleus) but also for protons and neutrons 
(using individual density distributions), which is listed in the 
Table II. For this calculation, we first find the range of the in- 
tegral, i.e., the lower and the upper limits of the axes (r and z) 
from the plots in Figs. 2 to 5, and then evaluate the integral 
for each case. 

Table II first shows that the number of protons and neutrons 
from the individual density distributins gives the total num- 
ber of nucleons which is calculated here from the total density 
distribution of the nucleus for each isotope of Ba with differ- 
ent fa solutions. Secondly, the internal cluster configuration 
of 112_122 Ba nuclei contain 12 C or 13 C as the most common 
cluster. For the oblate-state (/? 2 ~-0.35, Fig. 2), the density 
distribution is uniform for 114 ' 116 Ba isotopes, but in cases of 
H2,ii8,i20g a isotopes we find a 12 C or 13 C cluster, and for 
122 Ba a 23 Na cluster. The ground-state density distributions 
(around fa -0.25, Fig. 3) of 112 ' 114 Ba contain 13 C cluster, 
but 116_122 Ba contain clusters of some other elements like 
6 Li and 12 H. The density distributions for the prolate-state 
(fa —1.2, Fig. 4) are again uniform for the isotopes 112 ' 116 Ba, 
but show 8 Be cluster confuguration for 114 Ba. Note that there 
are no solutions for fa — 1.2 for the isotopic chain 118_122 Ba. 
For the very large fa ~10.5, the internal confugurations of 
the isotopic chain 112_122 Ba are of the form of two separated 
(identical) nuclei which are connected by a neck like confugu- 
ration. The neck confugurations contain simply the hydrogen 
isotopes 2 ' 3 H or 4 He nucleus. As already noted in the Intro- 
duction, the existence of 12 C cluster inside the Ba nuclei has 
been of interest both from experimental and theoretical points 
of views. The important point is that 12 C cluster is formed 
inside the Ba nuclei, and are not from the neck region where 
2,3 H or 4 He nuclei are shown to exist. In other words, 12 C 
constitutes the cluster structure of Ba nuclei. 



IV. SUMMARY AND CONCLUSIONS 

Concluding, we have calculated the gross nuclear proper- 
ties and the nucleon density distributions for the isotopic chain 
112-1 22 g a usm g deformed relativistic mean field (RMF) 
formalism with NL3 parameter set. The gross properties, like 
the binding energy, deformation parameter fa and the charge 
radius r c , show qualitative similarity between the experimen- 
tal and RMF calculated values. Analysing the nuclear density 
distributions, we get the internal or sub- structure of clusters 
in Ba isotopes which we find to consist mostly of 12 ' 13 C and 
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several other clusters like 12 H, 6 Li and 8 Be. Some heavier 
clusters of Na, CI, Mg, Ar and Ca are also obtained. With the 
increase of deformation (/?2 ~ 10 or 1 1), the B a nucleus breaks 
(fissions) in to two symmetrical fragments, releasing from the 
neck region some clusters of hydrogen isotopes 2 ,3 H or 4 He. 
This is an interesting result of the RMF(NL3) technique for 
nuclear structure physics. 

Clustering is also important for the decay of excited com- 
pound nuclei formed in nuclear reactions, where nuclei in the 
neck region could have correspondance with the measured 
fusion-evaporation residues consisting of light particles with 
Z <2. We know from the Introduction above that todate the 
decay of 116 > 118 > 122 Ba* compound nuclei in to intermediate 



mass fragments (IMFs), and symmetric and near- symmetric 
fission fragments are measured HH E2, EH EH EE[], but the 
fusion-evaporation cross-sections are not yet measured. 
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